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Abstract Wilh a standard deviation ot large at 2 em, the sea 
sate is (SSB) has become the dominant suce of error it 
Satelite altimetry. The operational SSB correction modes are 
{Werdimeaioaal 2D) empirical (paramelre or nonparametric) 
todels haned on Ue alneterensured wind speed (U) and 
Suncat wave height (SWH). However, these 2D SSB models 
‘hn mot entirely parameterize the range has variability, The SM 
Uiceraity may be lowered through improved SSB models 
Including aditional measurable or predictable coctis This 
paper proent a method fo climate the SSB from crsver 
flens by wing  theeditindona (D) timparsmetic 
"model. The modal i based om U, SWH fom Jason? altimeter 
acean observans, and the mean wave period (MNP) Irom the 
Earopean Centre for Median Bange Weather Forecasts 
(EEMWE) reanalysis project EX terim (The SS model 
eveloped withthe method presented in this paper b cale "+D 
SSB mde andthe SSB estimated with (be 2D SSB model i 
led D SSB cima Simulations indicate thatthe Wave 
period can greatly affect the SSB. Evaluated by the separate 
ual data Sets frm 2009 lo 201, the D SSB estimatet car 
increase the explained variance by 1-32 en, or LIS cm RMS 
‘lative othe traditional 2D SSB estimates based on U and SWH, 
Spatial evaluation of improvement shows thal the YD SSB 
‘lias are better than the traditional 2.D SSB estimates t all 
Tatts The enhancement frm 2D ta 3-D SSB estimates af 
‘eat dgulicane to improve the precion of (he alimeler 
product 


Indes Tersan state bas, rador alineler, crossover 
erences, thee dieasinal nonparamerie model, mean Wave 
period 


L Istropucrio 
za level sise, mainly caused by global warming, has 
scd mere and more atemton. One of the main 

purposes of radar altineters is to measure the mean sea level, 

However, the sea level measured by radar altimeters i lower 

than the te sea level because wave troughs are beter 


reflectors of radar energy than wave crests [1]. This effect is 
known as the electromagnetic bias (EM bias). There are two 
‘other sectae related biases: the skewness and tracker biases 
[2]. The EM bias, the skewness bias and the racker bias are 
generally called the sea state bias (SSB). Due to improvements 
im the precise orbit determination. technology and other 
seophysical corrections, the sea state bias is now the dominant 
ror term in satelite altimetry (3). 

Ta 1971, Jackson first discovered the phenomenon ofthe EM 
bias Trom radar observations in a tower experiment (4). During 
the past Tew decades, various investigators have ted to ue 
"oh geometrical optics [49] and physical opties [10-12] 
scattering theory to describe the EM bas. Numerous tower and. 
airborme experiments have been conducted to determine the 
Behavior af the EM bias [1-17] Nevertheless, itis sill hard to 
completely understand the physical mechanism of the EM bias, 
30 the estimation ofthe SSB still largely relies on empirical 
models. 

Empirical SSB models are usually developed with a 
parametric [2 18, 19] or nonparametric [20-15] method. 
Operational SSB models are generally based on the sea surface 
height (SSH) differences between repeat altimeter tracks, either 
crossover [2] or collinear differences (18). More recently, the 
sea level anomalies (SLAS) obtained by subrang the mean 
Sea surface (MSS) from the instantaneous SSH measurements 
have been used to directly estimate the SSB [2226] Empirical 
SSB models are classically developed by using the 
altimeter measured wind speed (U) and significant wave eight 
(SWH). However, these two-parameter SSB models can not 
entirely parameterize the SSH variability atributed to regional 
complexities in ocean wave climate (27, 28]. Backsatering 
coefficient. pseudo wave ae. swell height. mean wave period, 
RMS long wave slope and inverse wave age Were used 10 
esimate he SSB I9, 23, 25]. Tran presented a 
three-dimensianal (3-D) stellte SSB correction model in 2010 
and the results indicated that the new model could reduce by 
Sound 7.5% the total altimeter range error budgets for Jason- 
and Jason-2 [24]. In addition to U and SWH, Tran used the 
mean wave period ( Tm ) estimated from a numerical 
wind-wave model, NOAA'S WaveWatch3. WaveWatch3 was 
‘un on a global 1°x1° grid at 6 hourly time step in Tan's study 
The mean wave period is defined as 

Tm=m m, o 
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[[suor'auo o 
tem, represents he sti moment derived m te 
directional wave elevation deny spectrum. SU/.D. with 
fraguncy fad wave poppin dicun d. 

“ran developed he 3D SSB model using he SLAs dala sets 
StAsaccbtine by subuacting the mean se surface iom he 
tmansecus SSH marone Tia an m nice i 
sencrly bined fom a mean sea surface model Therefore, 
"rea wo psit sores of eror Tor the SLAs dt te 
Sos inte SSH meisremcts andthe sie etri the 
mean sea sce model. Any munge or geophysical correction 
Corinth SSH messe at dct iet n the SS 
Soma, vira e SSH difleencesexination is mue 
ly t roce siena enor hy repe tack andy in. 
2004. Laboc. [22] compared. the SSB cimus wing 
tometer demo, «inet diras ual SLAB 
Teqetordy. Tie eror anf rli eres give he 
aerei provided atthe SSH arid ints way 
However, te SLAS estes are significantly dien fom 
fe collber or conover estates, Io aon, the SS 
le iv fem ils nen n sace dl ae 
quie ferc and ber no erit EERO fr mc 
roce Thef, tne crossover nd colinear estimates are 
prefered forthe moment, but he SLAS eimai ni 
Someta ye 

Inthe paper. we present a method to estimate the SSB from 
cromover difrences by wing a theediensinal Q.D) 
onpasnori model Lie what Tran did in 2010, we also oe 
te Win speed te ignit wave height nd ic mean we. 
eod athe thee inp vaidles The ECM global 
omis product ERA-Imerim provides the mean wane 
rod ney The mean wave perd from ERA tm is 
Based on -1 petal momen (v Fn) 29, which dierent 
fom ht in Tras model, However, he 3D SSB model 
‘ing the mean wane evo ee provided by ERA nem. 
o has god performances, the 3-D SSB model will be mach 
re pai Beams the ECMWE rende product hae 
cen widely wed inthe alimder data pocceing an it 
‘alte over a ong me period 

"The sest ofthe paper i organized as fellows- To analyze the 
iei henya te sure peod and SSB, we fia 
‘Soins profs of th sen mc a dep wate vh e 
foarthorder Saks wave theory in Section Ht Section MI, we 
‘describe he onpurnetrc ciao method for the SB 
nimata. Estimating the SSB fom he SSH teens by 
‘sing 4D nonparametric model pas not been implements 
Score Comput with the 2D SSB estimation, te 2. SSB 
iun rogi mech mere companion tine, We pesa. 
the methods to reduce the computation tine in dis ation 
Sei 1V ines he dua sce andthe daa collation 
proce tn Sesion V, ye otras te SSB of Jn 2 ag 
fo 2D and a D nonparnati models. Zen-adjutment 
Brocesing forthe SSB mel nln the mean wae eid. 
cues i th scion, Son V k the conlason of the 
paper 


JL THERELANONSHIP BETWEEN THE WAVE PERIOD AND THE 
SSB 

Traditional operational SSB models are based on the wind 
speed and the significant wave height, Trans 3-D SSB model 
presented in 2010 indicated thatthe variance reduction of 
residual sea surface height (RSSH) was inreased ut all 
Jatiudes afier be mean wave period was introduced 10 the SSB. 
models. However, there have been no publications that describe 
the relationship berween the wave period and the SSB. In this 
section. we aim to discuss the relationship between the wave 
period and the SSB. EM bias in the SSB is mainly caused by 
nonlinear sca surface. Therefore, any factor elated o the 
nonlinear sea surface can affect the EM bias. 

‘Sea wave i generally studied from two aspects, On the one 
hand. uid mechanics i» used to study the movement of 
particles in the liquid and the two main theories are the Linear 
and nonlinear sea wave theories. On the oiher hand, sea wave is 
Considered as a random process and the movement of particles 
in the liquid is described statistically. The mean wave period 
used in this study is actually obtained in the second way. 
However, our main objective is to discuss the relationship 
between the wave period and he SSB, so we use the linear and 
nonlinear sea wave theories în this study. 

According to the linear and nonlinear sea wave theories, sea 
wave is linear when the waye heights small enough compared 
With the wavelength. However, linear wave theory is not 
suitable when be wave height is no longer small enough. 
Nonlinear sea wave theory mainly consists of the Stokes wave 
theory and the Cnoidal wave theory. Fig. I presents the ranges 
of validity for diferent wave theories [301 The Stokes wave 
theory is applied to deep water, but in shallow water, the 
Cnoidal wave theory s preferred. 


Most valid measurements of alimeters are in desp water 
Where the Stokes wave theory is preferred. For the Stokes wave 
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theory, he higher the order is, the wider the range of validi is. 
In this study, we use the fourth-order Stokes wave theory. 
According 1 the fourth-order Stokes wave theory, Ihe wave 
profile can be writen as) [31]. 


Duties 


st sa bee vns 


3 Ls. 
23a conatays Let cot 
i +L co, 


o 
ra o 
E 9 
[2 i) 
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where LT, H, k are wavelength, wave period, wave height and 
‘wavenumber respectively. The term a is an undetermined 
constant associated with the wave height If we know H and T. 
‘we can obtain the wave profile using expression (3}(6) 


AAKARS 


dir, a The vane hag ed we tbe wan eed anal) he 


Fig. 2 presets the wave profiles of the sca surface that are 
obtained withthe fourth Stokes wave theory. When the wave 
height is fixed. nonlinearity is weakened with the increase of 


the wave period On the contrary, when the wave period is fixed. 
nonlinearity is weakened with the decrease of the wave height. 
The wave profile of the nonlinear wave is usually an 
asymmetric curve with seep cress and flat woughs. The 
strength of the nonlinearity is mainly determined by the relative 
wave height H/A in shallow water, where d is water depth. 
However. in deep water, the stength of the nonlinearity is 
mainly determined by wave steepness & = H 7L. The larger the 
wave stepness is, the stonger the nonlinearity ix The 
waveleagth generally increases with the wave period. 
Therefore, similar to the wave height, the wave period can 
srcaly affect the SSB. When the wave height is xd, if the 
‘wave period increases, the steepness willbe decrease, 80 the 
magnitude of the SSB will decrease accordingly 
The wave profile g(x) in (3) satisfies 


["eo-dae tiet. o 


‘Therefore, the level z, = Lug +K a" forms the free surface on 


3 
Which the water particle is still, that is to say. the origin of 


cots s, Lat «Ef si i it vo 


surface, 
"The height of the wave crest with respect to the origin o£ 
coordinates is 


TM 
uy 


laesae ds 


The height of the wave tough with respect t the origin of 
coordinates is 


eu 
* 


Hr d wf catt ea a e 
E eer 
Prod o 

See deu. 

ants 


The half wave height with respect 10 the origin of coordinates 


i 


H,- kins) 


i AEN I. id 
Lie fw et! eats av 
3 BC 
‘Therefore, the half wave height with respec to the still water 
surface 


H =H- 


an 


Vue My s 
DuC. 

H, can be seen as the theoretical EM bias derived from the 
Stokes theory. We can obtain H, from (4). (5) (6) and (11). 
Fig. 3 shows H, that varies ong with the change ofthe wave 


height and the wave period respectively. When the wave period 
3s fixed, H, isan increasing funcion of the wave height. On 


the contrary, H, is a decreasing function of the wave period 
When the wave height is fixed 
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TIL NONPARAMETRIC ESTIMATION METHOD FOR THE SEA 
SpacBus. 


Gaspar fist introduced the noaparametic (NP) method for 
the estimation of the sea state bias in 1998 [20] and then 
developed it in 2002 [2I]. From then on. nonpuramctic 
estimation method has been widely applied io the SSB 
cstimation. In this section, we discuss the nonparametric 
estimation method for the sca state Dias, Part A presens the 
general multidimensional nonpuramciie regression model 
‘Sing the local linear regression estimator. Part B presents 
ars D nonparametric SSB model developed with the 
SLAs data. Part C describes the basic principle of estimating 
the SSB from crossover differences by using a nonpatametric 
model. The differences between the SLAs estimates andthe 
crossover estimates are also discussed in this part. Compared 
with developing a 2-D SSB model, developing à 3-D SSB 
model requies much more computation time. Part D gives the 
methods to reduce the computation time. On the one hand. we 
adop a new iterative algorithm. for sparse least-squares 
Problems, On the other hand, we improve the lecal linear 
Tegression estimator described in Pat A. 


A. Muttisdimensionat Nonparametric Regression Model 
Assume that y âs a random scalar variable jointly distributed 
with a random vector = (x x.) The general regression 
model can be deseribed as 
yers, a» 
where ris a regression function and £ is a zero-mean ertor 
term. Nonparametric estimation theory provides various 
estimators for the regression funetion rix) -ELv |x]. 1E we 
have n observations (s), we can write the estimator F(x) of 
the regression function r(x) = EL y |x as 


Xs 


a» 


where (xx) is a weighting funcion that as to satisfy 


Yana 


Gaspars research has shown that the local linear regression 
CLLR) estimator is the est choice to minimize the bias on the 
‘estimate [21] Therefore, the local linear repression estimator is 
Widely used to develop SSB models. For any vector z in the 
neighborhood of x. a local linear model for the regression 
function 7 can be vriten as 


AER 


an 


ra) 


ET us) 


We can find the value of B= (Pq efi" that minimizes 


ae 


Bts-A- Sg Pk n 


The dene function K,(x-x) 


KTH*-x)] is the 


m 
function that describes the weighting scheme for the 2, (x-x) 


weight term given in (13). The kernel function K assigns 
weights as a function o a scaled distance H*(X-x between 


‘sand. [His he dteinat of, a (px pyntndagilar 
scaling ma, generally al bandwith nt. The eel 
funcion K a nonnegative scalar funcion normalized 10 
saty J Roo 
pechailty demiby func. In this paper we used de 
Spherical Epanechaikov kernel, which proves o minimis the 
em iuis emor of be LR ema (2) 


Te is usually chosen to be a symmetric 


Fora 2-D SSB model based on U and SWH, x= (L, SWH), the 
kernel reads [21,33] 
Kaisa) ^n 
= an 


ERU 

E MH e, 
p" 

For a RD SSB model based and U, SWH and MWP, 

(U, SWH, MVP) tbe kernel reads [33] 


A EP P 1a 

Eana R T 

ER 
ee ee 


where bh tg) ie bandwidth vector, We wed 
local bandwith inti paper. Dea can he found in the 
Ines [21] The bandh vector hay ee as 
Rha 8) 
jiener iri m9)? eas a» 
Tf: ons eJ S) 0. 
where (s Pai a reference Banh We sls 
(haha, value Cis 09m 2s referred to the 
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‘erature [21]. nts) is the number of measurements in the box 


where xis fond and Ris the mean value of ni). 


Equation (16) ix a typical weighted least squares problem and 
the solution is: 


NEWN, NDW e» 
where Wdias(ycoX-X)) is the (rem) matrix of the 


EL 
" EL 
weighs, y* obs. Xp =); x 
1-3 


The so-called local linear estimator of F(X) is f, while the 


other components of are estimators of the first order partial 
derivatives of r [32]. We can thus write 


huata) =p, TOGWA, V NEW, e 
Where e, is [pest unit vector: ef 21,0...0). Sening 
T-éADWAU NSW, the weighing function can be 

2 8-3) T0. e 


Obviously. Fu conforms to (13) and, with a ile algebra 
,(x-X) in (22) can be shown to satisfy (14). 


B. Tvan's Three-dimensional Nonparametrie SSB Model 

Trans tree-dimensional parametric SSB model is based an 
‘he wind sped, the significant wave height fom Jason-1 and 
the mean wave period estimated from WaveWatch3 [M] 
Extracting the SSB from the SSH measurements is the key 10 
develop a SSB model. Tran extracted the SSB by subtracting 
the mean sea surface height from the instantaneous sea surface 
height. 

"The sea surface height measurement uncorrected for the SSB 
is denoted as SSH . SSH contains the geoid (h, the dynamic 
topography (7) the SSB and some measurement noise (o). 

SSH =h, 474 SSB +00 e 

1E SE, is the estimated SSB from a SSB model, the sea 

surface height measurement corrected fr the SSB (SSH ) can 
be writen as 


SSH =h, ope SSB-SS e» 
At cach location, the mean sca surface can be obtained by 
averaging the multiyear measurements t be location k. 
p 
es 
"The sea level anomalies can be obtained by subtracting 
the mean sea surface height from the instantaneous sea 
surface height. 
Ah, = SSH MSS, 
hk eh hme os SSB, 


eo 
SS -SSB, >, en maro 


Assume that all terms except for (SSB), in (26) can be 
considered as a single noise term 4, of zero mean. Equation (26) 
can be simplified: 

Ak, =(S8B), e en 

Assume thit the SSB is a function g of x, a vector 
containing U, SWH and MWP. Noting y, — A, then 

Ak, = SSH, - MSS, = (SSB), e - ex en. CR) 

"This is the typical expression o a regression model. Given a 

Set of n measurements (y x). where vis the residual sea 


surface height and x, = (4.420%) = (U. SWH, MWA) is a 
Vector defined by the measurements of U, SWH and MWP, the 
estimated SSB p(x) can be obtained from (21). 


C, Estimating the SSB from Crossover Differences by Using a 
Three-dimensional Nonparametric Model 

Unlike Trans 3-D SSB model developed with the SLAs data 
the 3-D SSB model we developed in this paper is based on the 
crossover differences, 

AL crossover poins, the geoid signal can be eliminated by 
forming differences between two SSH measurements taken on 
the ascending are (index 1) and the descending are (index 2) 


respectively. 
SSH, SH, m 
Assume tiat) o. - can be considered as a single 


the SSH 


2ero-mean noise term. The observation equation for 
Aliferences can be writen as 

SSH,- SSH, = SSB, SSB, +e. a» 
Assume ihat the SSB is an unspecified function g of 
(U. SWH.MWP) : SSB= $x). noting y 
the observation equation (30) can be writen as 


Y= etx) elm) +e. on 
IF x, -xis given, the conditional expectation of y i 
Alyx: =x1= 900) —Alek) |= 02) 


With a set of observations (v, Xx 
conditional expectations in (32) using 13): 


Soss-x Ern G-x) 0m 


Fora givens, if the value of pix, ) is known, (33) provides a 


we can estimate the 


em) 


practical method wo estimate (x). To estimate it, we can write 
(G3) for x equal to any of, 


PECES 


(Orin mati form: 

[zr [I 
Where His a (a>) identity matis, A is a (nxen) mamie with 
m 


X, -x,) expressed by the local linear 
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estimator of equation (14), gj 


TI 


E 


m 
The matrix E-A is singular because the rank of T= Aso 
Asdescibedin (30), only SSH differences are observed, so the 
‘SSB can only be determined to within a constant. Therefore, we 
have o impose the value of one oF the components of, say, 
ix) 94. S) cun be rewritten as 

Bh =Ay-Bya. ao 
where! =[ 0%), (B B) is a partition of T-A 
with B, the fist column of I-A. We now havea system of n 


equations with n-1 unknowns. A linear least squares solution is 
porem 


9 BIB, ira en 

“hit provides the missing estimates of gix, 72... 
“These can then be plugged iio (33) to estimate gi) for any 
value of. The sonprameti SSB esmalon problem i tus 
formally solved 

"The main difeence between te SLAs estimates and the 
ceromover eer the vay to ec the SSB. The former 
tract the SSB from he sea level anomalies while de later 
dors om the SSH ferences he cromover. Compared 
with the SLAS SSB estimates, the cronaver SSB estia 
Bave some advanage: 

1) Sone sync errors can be eliminated by the SSH 
<aerences teats they are unchangeable for a sume postion 
While ch sae errs cant be eliminated byte SLAN 

2) The SLAS ae obtained by subacting te mean sea 
arce from the tntneots SSH mensrenens. The mean 
ea surface ia ally obtained ora mean se surface model 
Tua is chaise from several aimee, Same enun: muy 
Spear since dier alim: we dilrmat methode 10 
‘process e da In addition, diferent mean ses surface models 
1 ped quine different SSB estimates 22 


D. Methods Jor Reducing the Compton Tome 

As shown in (37). for each x it requires large number of 
matrix operations, so the nonparametric method can hardly be 
applied to the real time processing of the altimeter data directly. 
To solve this problem. we can make a lookup table and estimate 
the SSB at each grid in the lookup table. Compared with the 
2-D SSB lookup table, the 3-D SSB lookup table has many 
more grid points. Therefore, developing a 3-D SSB model 
equites much more computation time. In this section. we 
describe the methods we used to reduce the computation time 
for developing à -D SSB model. 

1) The application of LSMR: BEB, in GT) is a 
lor-I«(r-D] max. Imvening the max requires a 
prohibitive amount of computation time when nis very large- In 
Gaspurs study in 1998 [20]. 500 crossover poins were 
randomly extracted for each cycle, so the dala subset was 
init to 9 = 500. In 2002, Gaspar [21 sed the LSQR (Least 
Square QR decomposition method) [3] algorithm to solve 
uation (36) and greatly reduced computation ime, so all 


valid crossovers in each cycle were used to develop the SSB 
model. In our sudy. we replaced the LSQR algorithm in 
Gaspars research wih the LSMR (Least Square 
‘minimum residual method) [35] algorithm 

Both LSQR and LSMR are iterative algorithms for sparse 
lincar equation and sparse least squares. The two algorithms 
find a solution x to the following problems: 
Unymneric equations: minimize ff subject tn Ax: 


on 
Linear least squares (LS): minimize Ax-w. 29) 


H E 
Di 


Regubrized least [4-4 
á Ei E N 
where Aet" eL", and mnor m2. The matin A 
is normaly large and sparse. 
The praca! sopping rules involve thee dinensnlse 
quits ATOL BTOL and CONLIN 


SI: Stopi Je. 8TOLj+ ATOL JAI |; an 
S2: Stop re ATOL Jall I [I 
S3: Stopif eond(A) = CONLIM. am 


The fü wo rules are applied 1o compatible and 
incompatible systems respectively while the third rule work in 
both systems. 

LSQR is equivalent to the conjugate gradient (CG) method 


applied tothe normal equation (ATA eye Ab. It can 
reduce s monii. where r, zb- A, is the riui 
for te approximate solution x, LSM s 
equivalent tthe minimum-residual (MINRES) method applied 
to de somal equation, so An] are monooniciy 


decreasing. In practice, fi also decreases monotonically and 
is never very far behind the corresponding value for LSQR. 
Therefore, although LSQR and LSM ultimately converge to 
similar points, it is sufer to use LSMR n situations where the 
Solver must be terminated eariy. 

After using the Epanechnikov kemel (36) becomes a sparse 
linear system, for which both LSQR and LSMR can be used, 
B,. and Ay—Byo, are comesponding to A. x and b 
tespectively in (38) - (40), We used both LSQR and LSMR to 
solve (36). LSMR produces results abit fuster than LSQR does. 
‘Alter replacing LSQR with LSMR, the computation time can 
be redoced by 105-156 


2) Improved local linear repression estimator: The inl 
ea les regression estimator of 3) is 
fia) =, - d OGWX, NEWS «ee X, is a (xp) 


matrix and Wis a (nn) matrix. When mis large enough, his 
requires a Jot of computation time. It is not obvious for a 2-D 
SSB model, but fora 3-D SSB model, much more computation 
time is required because there are many more grid points in a 
3D SSB lookup table. We atempted o avoid high dimensional 
"matis operations and made some improvements on the 
traditional ocal lincar regression estimator. 
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Fora 3D SSB model based on U, SWH and MWP, 
(U. SWH, MW), for any vector z of thee variables, à 


local incar model for the regression function + can be written 


rüb- Aes x els oko Ao (4) 
We have o find the value of B= (fi ff" that 


z-Ss-p. 


1B.“ )- Rs Kta) 
a9 


Bis 


‘To get, we can calculate the partial derivatives of Z and make 


hee az E 
em equal t zen. E. 
[3 
Four equations are obtine 
Say Sal + Sm: + Saa, 
Sat af Sil Sa P 
Sond Sd Sd Su 
sehn Soft Sf 
bere bc d and Sate writen as 
am Sy culty 
b- $na- a-s 
e= Èske) a 
4S yaya )K nx 
sy n) 


[see Sin Sue Sar] 
We seca] 5m San Su Sm) (ae) can be written in 
|S Sin Sen Suu 
Se Sus 
^ 
meni am. Sp" |. ten 
a 


Moore-Penrose pseudo inverse matrix of matrix $. 
The Local linear estimator of r(x) i 


COACE SN 


a 


where eF =(L0,0.0) , setting Quest s 71x) can be 


7.0004 Q2)-b4QG)-€+04)d 


Meo oce fett 


CERI 
a» 
The weighting funcion can be writen as 
prm 
[OW QU), 22) [I 
QUI, 3) QU), A) 


As shown in (22) and (50). the form of the weighting 
function obtained with the improved local linear regression 
‘estimator i different from that obtained with the traditional 
Jocul linear regression estimator, For a 3-D nonparametic 
model, he weighing funcion is oblained from 
T=eFOGWR, J XSW cei a (4x1) unit vector, X, 
(n) maix and Wis à (nxn) matis, XSW is a em) 
mars and XWX, is a (4x4) matrix Therefore, de 


multiplication of Xj and W requires dn’ multiplication 


operations. The muliplicion or XZW and X, requis 4' 
malipicaon operations. The inversion of XEWX, requires 
4° melipiaion operations. The mulipication of 
ODWASU amd XIW gemis d'n multiplication 
operations. The multiplication o'e} and (XEWX PEW 
requires Am mahiplication operations. Thus, the dme 
Complexity of obuining the weighing faction wing the 
‘etna local linear regression estimator is n°) 

Ax described in (7). computing a. b.e, dand 5, requiren: 2n. 
2m 2a and (mekétiln multiplication operations 
reel. S is a dei) matie. Therefore, de tine 
complexity of obuining the weighing function using the 
improved Vaca linear regression esümabr i» Oi). In 
comparison with the tational lc liar regression 
‘Estimator, te improved focal incar regression estimator can 
ely reduce the computation tine 


IV. DATA SOURCES AND CoLLOCATED SET 


A ERA-Interim Reanalysis 

The mean wave period used in our study is obtained from the 
ERA-Imterim reanalysis. ERA-Interim is the thi global 
reanalysis project produced by the European Cene for 
Medium-Range Weather Forecasts (ECMWF). ECMWF hus 
spared a efforts to produce global reanalysis far the state ofthe 
aümesphere, land and ocean during the last few decades 
ERA-IS and ERA-40 ar the two major global reanalysis 
projects, ERA-Interim is the successor of ERA-15 and ERA-40. 
‘The objective of the ERA-Interim i to prepare Tor a new global 
teanalysis project t replace ERA-40 gradually. ERA-Interin. 
provides reanalysis data from 1979 to the present time and 
Updates in wear ral time, Contasted wiih ERAW1, 
FERA-Iierim has some improvements in data assimilation, 
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forecast model and observing system. ERA-Interim uses à 
spectral T255 horizontal resolution (T159 for ERA-40) and 
bins a reduced NI28 Gaussian grid with an approximately 
uniform spacing of 79 km (125 km for ERA-40). Therefare. 
ERA-Interim owns a higher horizontal resolution (T255, NI28 
approximately 071° ) han ERA-40 (T159, NSO approximately 
125"). With bilinear interpolation techniques the data i the 
NIS Gaussian grid is imerpolsted into different kinds af 
atiadefongitude grids from 0.125°t03". 

The ocean wave analysis in ERA-Interim incorporates an 
optimal interpolation method te constrain predided wave 
spectra using altimeter wave height observations. ERA-Interim 
ties reprocessed ERS-1 and ERS data Hom ESA, and 
near-real-time data from ENVISAT, JASON-1, and JASON-2. 
"The horizontal resolution of the wave model in ERA-Interim is 
110 ks and the wave spectra resolution is determined by 24 
directions and 30 frequencies [36]. ERA-Interim updates data 
Tour mes per day at 00. 06, 12 and IS UTC. A variety of data 
in uniform laiadelongitude grids can be obtained Rom 
ERA-Imerim datasets. The mean wave period we used is 
projected onthe grid oF 0.125" x 0.125 


B. Jason-2 Altimeter 
The wind spoed and the significant wave height used in this 
study are obtained fom the version "D' GDR of Jason-2. 
Jason-2 was launched in June 2008 and placed on the same 
ground tack as its predecessors TOPEX/POSEION (TP. 
launched in August 1992) and Jason- (launched in December 
2001). Jason-2 caries the Poseidon- radar altimeter and 
occupies non-su-synchronous orbit at an alitude of 1336km 
‘with an inclination of 66-15" Poseidon-3 is a dua-requency 
radar altimeter that operates at 13.575 Gl (Ku-band) and 5.3 
GHz (C-band). The repeat cycle of Tason-2 is about 9.91 days, 
which means that every location along me Jason? 
arouni-track is measured approximately every 9.91 days 


©. Collected Set 

Mean wave period is only avilable at 00, 06, 12 and 18 UTC 
every day, so we have to interpolate the daa onto the altimeter 
ground tack locations. In this study. bilinear interpolation is 
used in space and linear interpolation is used in time. 

"The crossover data we use in this study i from 2009 to 201 L. 
Gaspar computed the crossovers cycle per cycle, concluding 
that there are 6x10 to 7x10 valid crossovers per eyele. 
Unlike Gaspar [20-21], we computed the crossovers not only 
cycle per cycle but also between two adjacent cycles to get 
more crossovers. A group of crossovers consist of the 
crossovers from the curent cycle, the crossovers of the 
ascending arc for the current cycle and the descending are far 
the next cele, and the crossovers o the descending arc for the 
current cycle und the ascending are far the next cycle. Afr 
quality control we can get 15x10 to 18x10 valid erossovers 
Per group, which ar approximately three times that computed 
Cycle per eyele. There are few points that are completely 


consistent with the crossovers of the ascending are and the 
descending are. We select four points at cach side of the 
Crossover and use spline interpolation. (o interpolate He 
measurements 1o the crossover. 

“Gaspar computed the SSB estimates cycle per eyele and then 
computed the mean value of the individual SSB estimates. 
Unlike Gaspar. we fist put the crossover daa from all the 
groups together aad then randomly extract 8000 crossover 
points to compute the SSB estimates every time. Finally we 
Compute the mean value ofthe individua SSB estimates. In this 
Way, we may obtain more measurements from the data-ponr 
eas because there may be few measurements in the dat ot 
areas in a single group. 


V. ESTIMATION OF THEJASON-2 SSB 


A. Zero-adjusonent Processing for SSB Models Including 
wr 

As described in Section tthe rank of T-A in equation (35) 
Js equal to m-l and it may bring uncertainty problems. To 
iint the uncertainty, we have to impose the value af one of 
the components gy For a traditional 2-D nonparametric SSB 
model based on U and SWH, afier obtaining the estimated SSB 
eis) using (33), we should subtract oU =0.SWH =0) from 
XU.SWH) arall grid points 1o ensure thar the SSB is zero over 
a flat sea with no wind. We call this process “zero-adjustment 
processing’. Afer MWP is introduced to à SSB model, 
Zero-adjustment processing is very important. In this part. we 
discuss how to deal with zero-adjustment processing fr both 
2D and 3-D SSB models- 

Fora D SSB model based on SWH and MWP, we 


introduce two methods to do deal with this problem. In the finst 
‘method (2D method 1), we ensure that the SSB will be zero if 
SWH is zemo. H is implemetied by subtracting 
SW =O. MWR) Hom SWH. MVP) at all grid points 
The SSB is given in the form f 4 2-D look-up table, so at each 
‘ow, the obtained SSB ox) with (33) wil subtract the fist SSB 
to ensure XSWH =0. MWR) -0. In the second method (2D 
method 2), the obtained SSB estimates at all grid points subtract 
Same value, Fig presents the distribution histograms of U, 
SWH and MWP. The range of MWP in our 2-D look-up table is 
0-18 s and the mean of MWP in our data set is about s. Three 
values e(SWH =O.MWP=0) . (SWH =0.MWP=95) and 


SWH =0.MWP = I) are provided in this study. There are 
few measurements when MWP=0 or MWP=išs, so 
ASWH =0,MWP=9s) is more accurate compared with 
SWH -OMWP-O) and (SWH =O. MWP= IN) 


Therefore, the value @( SWE 
selected in this study. 


AMWP =95) is the value 
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Fora 3D SSB model based on U, SWH and MWP, 
“tero-adjustment processing is more complex. The SSB is given 
282 3-D look-up table. SWH, U and MWP are corresponding to 
‘he ro, column and page. Like what we do for the 2D SSB 
model based on MWP and SWH, we also introduce two 
methods o deal with this problem. For the first method (-D 
method 1), we ensure thut the SSB is zero when both U and 
SWH are zero. At each page, the obtained SSB fom (33) 
‘will subtract he SSB value that es in he first row and the first 
column to make o(SWH -0,MWP =0) =0. For the second 
method GD method 2), same value is subtracted from the 
obtained SSB at all grid points. Similar 1o the above 2-D SSB 
medel there are tee such values to. select fome 
0,SWH =OMWP=0) , gU =0.SWH » 
MW =I8s) . Similarly, there are few 

ss the obtained 


MwP 18) 
MIWP =9 is the value we select. 

We used the explained variance to evaluate different SSB 
‘models. Explained variance is the reduction in the variance of 
SSH differences obuined when applying the given SSB 
‘correction und it is a main criterion to assess the quality of 2 
‘SSH model [20]. Similar to what Tran did [30] we computed 
the explained variance using several complete-year data sets. 
"The year 201 version solutions are evaluated by not only 2011 
measurements but also 2009 and 2010 altimeter data. Results 
from 2009 version of the models are alo provided for 
comparisons 

For the 2-D (MWP and SWH) SSB models, we first 
calculated the explained variance using the global data. TABLE 
T presents the explained variance calculated with the data from 
2009 o 2011, The 2-D SSB model with method 2 has larger 
explained variance im all cases. We thea binned the data 
between 655 and 65°N from 2011 into 10° latitude bands and 
calculated the explained variance in each band, Fig. 5 presents 
latitude distribution ofthe variance explained by the 2-D SSB 
models with method 1 and method 2 respectively versus the 
D benchmark. The variance explained by the SSB model 


developed with method 2 is larger than that with method 1 at 
most Jatitudes. Therefore, the SSB model developed with 
method 2 is beter and we will use the second method in 
develop the 2-D SSB model based on MWP and SWH. 


Eas ETT 


opin trà Ib LTEM Geen ae Stal by 


For the 3-D (U, SWH and MWP) SSB models, we used the 
same way as the 2-D SSB models o evaluate the wo methods. 
TABLE I presents the explained variance calculated with the 
data fom 2009 o 2011 Similar to the 2-D SSB models based 
on MWP and SWH, the 3-D SSB model with method 2 is bener 
ata glabal scale. Fig. 6 presents the Iatitude distribution of the 
Variance explained by the 3-D SSB models with method 1 and 
method 2 respectively versus the 1-D benchmark. The SSB 
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model with method 2 is better than that with method 1 at all 
atitudes. Therefore. for both 2-D and 3D SSB models 
including MWP, subtracting a same value is a beter choice, 
ASWH =O MWP=95) and piU =0,SWH =0,MWP =93) 
are the values thar should be suberacted for a 2-D and +D SSB 
models respectively 


aE 


e A AET 

Preset con t AUE RE 
EPI An 
MP LO ue m ene 
cen a a pecs 


She hn ge SSA ag ee we 


‘To evaluate 2D US and 2D MS, we fist calculated the 
explained variance using the global data from 2009 to 2011, 
‘The SSB estimates are obtained with the erossover data whose 
time differences are less than 3 days. TABLE II presents the 
explained variance calculated with the crossover data of which 
the time differences are Iess than 3 days. The SSB correction 
should be added io the radar altimeter measurements of the sea 
Surface height (SSH). TABLE II also gives the variance of the 
Crossover SSH differences corrected for SSB. Ideally, the 
rossover SSH differences should be zero. However, because 
of ocean varisbility and some other conection errors, the 
crossover SSH differences are hardly zero. A better correction 
should decrease the variance of the crossover SSH differences. 
In fact, for two different models, the increase in the explained 
Variance equals to the decrease in the variance of the crossover 
SSH differences, As shown in TABLE IIl, compared with 
2D. US. 2D. MS inereases the explained variance by 0.38-0.59 
mf. The average of the six vales D. MS -2D. US) is 047 
em, or 0.69 cm RMS. The repeat cycle fr Jason-2 is about 
991 days. TABLE IV gives the results oblained by the 
‘crossover data whose time differences are less than 9.91 days 
Compared with 2D US, 2D MS increases the explained 
variance by 0.10-0.38 m The average of the six values is 0.21 
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cm. or 045 cm RMS. Therefore, the 2D. MS is beter on a 


slat scale. 


Therefore, MWP contributes more to the SSB at middle and 
: Jow latiudes iile at high Latitudes U contributes more. 
E] 


ndi er mcs edn 1a nds Rea 


To assess the performances of 2D_US and 2D MS at 
ferent regions, we binned the crossover data of which the 
time differences are les than 991 days into 10 Latitude bands 
and computed the explained variance in cach band. Fig 
presets the explained variance in comparison to à 1D 
(SASSWH) SSB model 2D. MS is bener than 2D_US at most 
Latudes while 2D_US is better only above S0°S and SDN. At 
middle and low latitudes, 2D.MS is obviously. beter 


C. SSB Estimated from the Crossover Differences by Using A 
ED Nonparamenic Model 

The 3-D SSB model based on U, SWH and MWP GD, USM) 
Js finally presented as a 3-D lookup table and the ranges of U, 
SWH and MWP are respectively 0-30m/s -12m and 0-135. U, 
SWH and MWP are all divided into S0 bands. "Therefore, the 
3-D SSB lookup table has SOSO SOgrid poins, As shown 
3n Fig. 4 the mean values of U, SWH and MWP are 773 ms, 
234 m and D14 s respectively. There are neatly no 
observations when any of the thee variables s beyond the 
range we give, so the ranges of U, SWH and MWP for the 3-D 
ookup ble are appropriate. In order to ter represent the 3-D 
SSB model, we present the correction values in the form of +D 
arrays by firing a ind parameter. Contours af the SSB values 
are displayed in em and shown in Fig. 9. The SSB lookup table 
3s obtained with one year of dta from year 2011. The fixed 
Values are near the mean values to ensure that there are enough 
‘measurements In our study. the three fixed value we choos are 
MWP= 95. U=75 mand SWH = 2 S respectively. Fig. (a) 
and (b) agree closely with Fig. 7 (a) and (b). When U and MW 
fe fied: the magnitude of the SSB is an increasing function of 
SWH. In contrast, the magnitude of the SSB is a decreasing 
function of MWP when U and SWH are fixed. In the data-rich 
region. there is an approximately 4 cm variation versus MWP. 
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ith bias decreasing as the wave period increases. As shown in 
Fig. 9 (c). when SWH is fixed, variations with MWP are more 
obvious than that with U in data-rich region. For a given U, the 
magnitude ofthe SSB is a decreasing function of MWP. 

As analyzed in Section IL in deep water the nonlinearity of 
‘he waye is mainly determine by the wave steepness, whichis 
the rai of the wave height and the wave length, EM bias, 
mainly caused by the nonlinearity of the waye, isthe main 
component of the sea state bias. EM bias increases with the 
Strengthening of the wave's nonlinear. Therefore, EM bias 
varies with wave steepness. When the wave height is fixed, the 
Seepuess is mainly determined by the wavelength. The 
Wavelength generally varies with the wave period linearly. 
Therefore, as long as the wave height is fixed, the wave period 
can determine the steepness and then affect the EM bias. With 
the increase of the wave period. the nonlinearity weakens. so 
the EM bias also decreases. Therefore, the resus in Fig. 9 are 
appropriate. 

Fig. 10 compares the theoretical $SB derived from the 
Stokes theory with the SSB derved from the empirical SSB 
models. For the SSB derived from the Stoker theory, the 
ascissas of Fig. 10 (a) and Fig. 10 (b are the wave height (H) 
ånd wave period (T) respectively. However, for the SSB 
derived from the empirical SSB models, the abscissa of Fig. 10 
{and (b) are respectively the significant wave height (SWH) 
and mean wave period (MWP) The variation trend ofthe SSB 
derived fom the empirical SSB models is consistent with that 
derived om the Stokes theory in the whole. When other 
conditions are be same, the magnitude of the SSB is an 
increasing function of the waye height and is a decreasing 
function of the wave period. However. the theoretical SSB 
derived from the Stokes theory is a Tew center lower than 
the empirical SSB correction im most regios. The lines 
determined by the empirical SSB correction with 3D_USM or 
2D. MS are more close to straight lines than that determined by 
the theoretical SSB conection implied by the Stokes theory 
oin Fig. 10 (aand Fig. 10 @). The mechanism of the SSB is 
uite complex, the theoretical SSB derived Rom the Stokes 
"hcory in this paper is simplified. Many Factors can affect the 
‘SSB and we can hardly describe the SSB only using the wave 
height and the wave pti, 

Since we have selected XU =0,SWH =0.MWP - 95) vo 
perform the zero-adjustment, there may be some positive SSB 
correction values. Fig. 11 presents the histogram of the 3-D 
SSB correction values. There are a small number of positive 
Values, the percentage of which is only 0.016%. The SSB 
should be negative in theory. However, only afer all other 
fers are well corrected can the SSB model be developed 
Some residual rors may affect the estimated SSB values, so 2 
small number of positive estimated SSB values are acceptable. 


iix arger than Ii amples 
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As with the 2-D SSB models, we use the explained variance 
to assess the performance of the 3-D SSB model. The main 
objective is to determine whether the performance of 3D_USM 
s improved relative to 2D_US and 2D. MS. We fist calculated 
the explained variance using the global crossover dta from 
2008 w 2011. For comparison, we also calculated the explained 
Variance using the SSB conection values in Jason-2 GDR. 
TABLE V presents the results calculated wit the crossover 
la whose time differences are less than 3 days. The change 


s 


trom 2D_US t0 3D_USM provides 1.25-1.40 em improvement 
im the explained Variance. The average of the nix Values 
OD _USME2D US) is 132 cn or LIS em in RMS sea level 
‘estimate improvement. Compared with 2D. MS and the SSB in 
Sinco-2 GDR product, 3D. USM can increase he explained 
‘alan by respectively Den and 1.61 em or 0.92 em and 
1137 em RMS. TABLE VI gives te results allied wih the 
crossover data of which the time diirence are less than 9.91 
days. Similar to TABLE V. JD USM provides the largest 
explained variance among ali the SSB models in TABLE VI. 
Compared with 20_US,3D_USM can decrease the variance or 
"he SSH ditterences by 1.10132 cn and the average ofthe six 
"ales QD. USN-2D. US) is 120k or 1:09 em RMS, The 
Average of the explsined variances increased by 3D-USM 
‘lative to 2D. MS and the SSB in Jason-2 GDR product are 
‘especie 099 en and 1.38 em”, or 099 cm and LIB em 
RMS. Therefore, the 3-D SSB modelis benter han the 2-D SSB 
models on a global sesle 
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To evaluate the performance of 3D USM at different 
latitudes, we binned the crossover data whose time differences 
are less than 9-91 days inno 10° atitude bands and computed the 
explained variance in each band. The explained variance 
Telatiye to the ID (3.8565WH) SSB model is shown in Fig. 12. 
Fig. 13 presents the mean values and the standard deviations 
(Std) of SWH, U and MWP respectively. Like Gaspar s results 
in 1994 [2 the standard deviations of both SWH and U 
increase with the Latitude, However, MWP differs a lot from U 
and SWH. The mean value of MWP is larger in the southem 
hemisphere while the standard deviation is larger in the 
orbem hemisphere. Both the mean value and the standard 
deviation of MWP are the largest at middle latitudes both in the 
southern and northern hemisphere. As described by Gaspar and 
Tran [2, 25] the variances explained by all SSB models 
increase towards the poles because the standard deviations of 
bonki U and SWH increase towards the poles, The explained 
variance at the equatorial is very small because of the small 
Standard deviation of SWH. The explained variance is very 
Bigh ar high atitudes because the standard deviations of U and 
SWH are very large. 

‘As shown in Fig. 12, 3D_USM is the best oT all the SSB 
models. Compared with 2D.US. 3D_USM is better at all 
atitudes. At most latitudes, the improvement is more than 1 
cn. The improvement is obvious at low and middle latitudes. 
Compared with the 2D MS, 3D USM also increases the 
explained variances at ll Tatitudes, especially at high latitudes. 
"This indicates that U obtains some SSB information that cannot 
be explained only by SWH and MWP. The wind speed can 
affect the roughness of the sea surface and then affect the radar 
ros section. According 10 the theoretical explanation of the 
EM bias, the general definition of the EM bias is the normalized 
correlation between the radar cross section and the long wave 
ivan [8, 9]. As explained in Miller's research [III 
non-Gaussian long-wave statistics and the hydrodynamic 


modulation of short-wave amplitudes are the two major 
contributors to the EM bias- The presence of wind- generated 
Small-wave roughness is crucial to the portion of the bias 
caused by hydrodynamic modulation. so the EM bias is 
cct by the local wind speed. As showa in Fig. 13. at high 
latitudes, the standard deviation of U is very high wile the 
standard deviation of MWP is relatively low. At these Latitudes, 
U may contain more SSB information than that MWP contains. 
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In onder o assess the performance of 3D. USM at different 
time intervals, we binned the crossover data from 2011 ino 
[day bands and compute the explained variance in each band 
The explained variance in comparison to the ID (-3.88%SWH) 
‘SSB model is shown in Fig. 14. The variances explained by all 
the SSB models increase with the me Intervals. The dynamic 
topography and tide varies with the time. When the time 
interval is large, the uctustons of the dynamic topography 
and tides wil be great therefore, the SSH differences vary 
‘really. As shown in Fig. the performance of 3D, USM is 
the best in all bands. The improvement from 2D_US 10 
3D. USM is more than | cm in most bands, The improvement 
is almost the same in all bands. There is a gres jump in 
explained variance going from 05 day to 1.5 day. There may 
e reo reasons, On the one hand, a the same crossover point 
within a day, SWH, U and MWP measured from the ascending 
pass are almost the same as that measured from the descending 
Pass. Therefore, the SSB difference is small at the same point 
nd the SSB differences make litle contribution to the SSH 
Aferences. On the other hand, ocean variability is very Weak 
‘within a day, so the variance of the SSH differences is very 
Smal. 


As discussed above, 3D_USM is beter than the 2-D SSB 
models at all lites and all time intervals. This improvement 
an be seen as a great improvement in comparison to 2D. US 
and 20 MS. 


D. Comparison been the Jason-2 Derived SLA with sit 
Measurements 

Seu evel anomaly (SLA) is an important parameter used in 
‘oceanographic investigations. SLA is obtained by subtracting 
the mean sea surface (MSS) from the insantancous SSH 
‘measurements, The MSS model CLS2011 is sed inthis ty. 
In order o know whether the 2-D SSB estimates can improve 
the accuracy of SLA, we compared the satellite derived SLA 
with some tide gauge measurements from Senetosa in Corsica 
Inland calirtionhalidtion ste In Senetsa, main calibration 
Site for TIP and Jason missions is located at 40 km from Ajaccio. 
Fig. 15 presents the general configuration of Corsica absolute 
calibration sie. There are three tide gauges (M3, MA and MS} 
submerged on both sides of pass 085 af Taon-2. M4 and MS 
re installed very close together (<30 em). in order to monitor 
data benter, 


We used the data from cycle 50 to eyele 150 af Jasom-2 
GDR. Fig. 16 presents the MS tide gange derived SLA, Jon 
derived SLA with 2-D SSB estimates and 3-D SSB estimates 
respectively. Jasoa- derived SLAs with both 2D SSB 
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estimates and 3-D SSB estimates agre closely with tbe MS tide 
[ge derived SLA inthe whole, TABLE VII presents Jason-2 
SLA bias fom tide gauges. As shown in TABLE VII, the 
variance of Jason-2 SLA bias from the tide gauges can be 
reduced by 0624,67 cm when replacing the 2D SSB 
‘mates with the 3-D SSB estimates 


E. The Performance of the 3-D SSB Estimates inthe Coastal 
Regions 

The open-ocean altimetry i mature and is a great success far 
satellite-based Earth Observation. However, the coastal zone is 
m important domain where altimeuy remains grossly 
underexploited. Changing cireulation, sea level and sea stale 
have by far the largest impact on human society. Revacking 
algorithms and the geophysical corrections are the crucial 
territories for coastal altimetry. Therefore, the improvement of 
the SSB corrections is very important foe coastal altimetry. In 
this part, we aim to assess the performance of the 3-D SSB 
‘estimates in the coastal regions. 


— 


Crossover SSH differences are he main ool to analyze the 
overall altimeter system performance. They allow us to analyze 
the SSH consistency between ascending and descending passes 


137, 38). We binned the daa between the coast and 100km 
offshore ino Skm-bands and computed the crossover SSH 
Variance differences (the variance of the crossover SSH 
differences with the 3-D SSB estimates minus the variance of 
the crossover SSH differences with the 2-D SSB estimates) in 
‘ach band. The results in the 0-10 km band are discarded due o 
the presence of land in the altimeter footprint. Results are 
shown in Fig. 17. The variance of the crossover SSH differences 
can be reduced in the whole coastal regions by using the 3-D 
SSB estimates In the 10-20 km band, the eduction of the 
variance of the crossover SSH differences fom the 2-D SSB 
estimates to the 3D SSB estimates is extremely high and it may 
be due To he lite number of the crossover points in these 
bands 


H 


H 


Ta the coastal regions, the number of the valid crossover 
points is relatively smali and the results may be not convincing 
from stia point of view. A mor robust assessment is to 
compare the along track SSH differences fr cach SSB solution. 
Every location along the Jason-2 ground track is measured very. 
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approximately 9.91 days. We binned the data between the coast 
fant 100km offshore into Sku-bands and computed the along 
"rack SSH variance differences (the variance of the along tack 
SSH differences with the 3-D SSB estimates minus the variance. 
ofthe slong tack SSH differences with the 2-D SSB estimates) 
in each band. The results in he 0-10 km band are also discarded. 
Resuhs are shown in Fig-18. Compared withthe 2-D SSB 
estimates, the 3-D SSB estimates can reduce the variance ofthe 
long tack SSH differences in the whole coastal regions. The 
‘eduction of the variance of the along tack SSH differences 
from the 2-D SSB estimates to the 3D SSB estimates increases 
with the distance rom the coast. To compare the performances 
Of the YD SSB estimates in The coastal regions and in the 
opea-ocean regions respectively we binned the daa between 
‘he east and 1000km offshore into 20kmebands and computed 
the along wack SSH variance differences in each band. Resulis 
ae shown in Fig.19. Similar to the results in Fig 18, the 
reduction of the variance ofthe crossover SSH differences from 
the 2-D SSB estimates tø the 3D SSB estimates increases with 
the distance from the coast, Because all the other conditions are 
the same, the improvement of the SSH is due to the 
improvement of the SSB corrections. 

As analyzed above, both the crossover SSH differences 
variance and the along track SSH differences variance can be 
seduce in the whole coastal regions by replacing the 2-D SSB 
‘estimates with the 3-D SSB estimates, This means that the 3-D 
SSB estimates can improve the whole allnetry system 

the improvement 


p 
improvement in the apen-occan regions, In shallow waters, 
‘water depth can greatly affect the nonlinearity of the sea surface 
and therefore affect the SSB. In the future, we will free study 
‘he SSB in the coastal regions and introduce he bathymetry or 
some other parameters to the SSB model for coastal ahimerry. 


regions is relatively small compared with the 


E. Comparison with the 3D SSB Models Using the Wave 
Period from WaveWatehs 

1) Comparison wit the -D SSB model using the wave mean 
period from Wave Watch The mean wave period Tran used is 
cstimated fom a numerical ocean wave model, NOAA'S 
WaseWacha (NWW3). NWWS was un on a global «17 
ar. a 6 hourly time step and ECMWF wind fields were sed 
To force the wave model in Trans study. The mean wave period 
we used im this paper is provided by ECMWF reanalysis 
ERA-Inierim dicey and we do aot have to run the wave 
model by ourselves The mean wave period (MWP) from 
ERA-Interim is based on -I spectral moment m., /). On the 
contrary, the mean wave period (Tm) estimated fom NWW3 in 
Tran's -D SSB modelis based on 1 spectral moment (a, 
‘Moreover. the mean wave period rom ERA-Interim is modeled 
withthe Wave Model (WAM). Both WAM and WaveWatch3 
ae the third generation wave models. Wave Watcha i in the 
Sprit of the WAM model but differs from WAM in in many 
important points such as the governing equations, the model 


secure, de 
paramcterizations. 

Tn the NWW3 products provided by the National Oceanic 
and Atmospheric AdminisatiowNational Center or 
Environment’ PreictiwNational Weather — Senice 
(NOAAINCEPINWS the mean wave period is not provided 
over a long time period. Therefore, if we want to use the mean 
Wave period to develop the 3-D SSB model, we have to run the 
WaveWateh3 model by ourselves. In this paper, WaveWatch3 
as run on a 0.25°x0.25" grid and at a one hourly time step. 
The mean wave period (Tm) we used is obuimed from 
September 2, 2014 to October 21, 2014 over the integration 
domain from TN to 45°N in Latitude and I0S-E 10 150E in 
Jong, Fig. 20 presents the distribution of Jason-2 ground 
tracks we used to develop the SSB models. There are only 2 
small amount of crossover points, so we used the collinear data 
Set to develop the SSB models in this part. The SSB models 
developed with the collinear datu and the crossover data are 
"denis from the methodology point of view. AI the 
‘measurements are interpolated onto a reference tack of Jason-2 
and the collinear differences are the differences between two 
open cycles The two 3-D SSB models using MWP from 
ERA-Imterim and Tm from WaveWateh3 respectively are 
developed with the same collinear data sei. Furthermore, o 
assess the performance of the 1Wo 3-D SSB models, we aso 
developed three diferent 2-D SSB models. The spatial 
resolution of MWP from ERA-Iaterim we used to form the 
collinear data is also0.25°% 02 

As with the SSB models develope with the crossover data, 
we also use the explained variance o assess the performance of 
dierent SSB models in this pan. Fig. 21 presents the 
corresponding explained variance by different SSB models 
tested across the collinear data. The 3-D SSB model using 
MWP from ERA-Interim displays the highest reduction with 
34.10 cni, The 2-D SSB model based on SWH from Jason-2 
GDR and MWP from ERA-Interim displays a lle higher 
reduction compared with the traditional 2-D SSB model based 
U and SWH from Jason-2. The 2-D SSB model based on SWH 
from Jason-2 GDR and Tm obtained from WaveWatch3 yield a 
decrease in toa variance of 32.63 emè, wbich is 0.18 em lower 
than that explained by the 2-D SSB model based on U and 
SWH from ason-2 GDR. In Trans stdy in 2006 [25] the SSB. 
model based on SWH from Jason-1 GDR and Tm obtained 
from WaveWatch3 also yields 0.4 emè explained variance 
over than that explained by the 2-D SSB model based on U 
and SWH from Jason-1 GDR. However, the 3-D SSB model 
used on U and SWH from Jason-2 GDR and Tm from 
WaveWachā displays higher reduction than the aditional 2D 
SSB model, Therefore, both the two 3-D SSB models using 
MWP from ERA-Imerim and Tm estimated from Wave Wateh3 
respectively are able to reduce the total variance in the SSH 
Aiferences relative to the traditional 2-D SSB models. 


numerical mede and the physical 
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Fig, 24. Redacon in vance of SEN ems perm by he ie 


2) Comparison with the 3-D SSH model using the peak wave 
period from WaveWatch’: The mean wave period it mot 
provided over a long time period in the NWW3 products from 
NOANINCEPINWS. However, the peak wave period (Tp) is 
directly provided in WaveWatch3 venion 222 hindeast 
‘reanalysis product, which provides the peak wave period from 
‘Aug 1999 up to nov. We sim to develop a 3-D SSB model 
bused on U and SWH from Jason-2 GDR data and Tp from 
WaveWach3 version 222 hindeast reanalysis product using 
the crossover data set. The peak wave period we used 10 
develop the 3-D SSB model is on x 6 hourly temporal and 
(0.25 spatial resolution which covers he globe from FS t0 
D 

TABLE VII presents the magnitude of the explained 
variances obtained with several SSB models. Fig. 22 presens 
Taide diibuion of the variance explained by several SSB 
models respectively in comparison io à ID (-385SWH) 
benchmark. As shown in TABLE VIII and Fig. 22, the 3-D SSB. 
‘model using the peak wave period from WaveWatchs version 
222 hindeast reanalysis product is better than the traditional 
D SSB model based on U and SWH. However, the 3-D SSB 
‘model using the peak wave period is not as good as the 3-D 
SSB model using the mean wave period from ECMWF 


teanslysis ERA-Inerim. There may be to resons for the 
dierences between the two 3-D SSB models. On the one han, 
the mean wave period fom ERA-Imterim is obtained from 
WAM model while the peak wave period from WaveWateh3 
Version 222 hindeast reanalysis produet is obtained from 
WaveWatch3 model. On the onher hand, the influence of the 
mean wave period on the SSB may differ from that of the peak 
Wave period on the SSB. 

As discussed above, the 3-D SSB model using the mean 
wave period or the peak wave period has positive impact on 
seducing the variance in sea level measurements relative 1o the 
ttadiional 2-D SSB model. The mean wave period can be 
defined in many ways and depend on diferent specual 


[I 


on 


where m, =[|SUOI AMO represents de statistical 
moments derived from the directional wave elevation density 
spectrum, SC), with frequency f and wave propagation 
direction Ü. The mean wave period Tran used io develop the 
2D SSB model is T,. The mean wave period from 
ERA-Imerim i7, ,. Future efforts will tempo develop 3D 


SSB models using different wave periods and find ou which 
‘wave period contibues most to the SSB. 


Fig 22 Late dian a e varane Cane by several SSB mee 
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Operational SSB models are based on altimeter-measured U 
and SWH. In this study, we estimate tbe SSB fm the 
crossover differences by using a 3-D nonparametric model 
based on U, SWH fom Jason-2 GDR data and MWP from 
ECMWF reanalysis project ERA-Interim. 

We first simulated the profiles of the sca surface in deep 
water with the fourh order Stokes wave theory to study the 
‘lationship between the wave period and the SSB. Results 
Show that the wave period may greatly affect the SSB. This 
provides theoretical support for introducing the mean wave 
period wo the empirical SSB correction models. 

Compared with a 2-D SSB model based on U and SWH, the 
lookup table far a 3-D SSB model based on U, SWH and MWP. 
has many more grid points. Thus, the 3-D SSB estimation 
requires much more computation time, We introduced LSMR 
iterative algorithm and improved the traditional Toca linear 
segesson estimator, so the computation time was greatly 
reduced. 

After inducing the mean wave period o the SSB models 
developed with nonparametric estimation method, we should 
Pay atieton o the 2ero-adjustment processing. We gave and 
compared two methods for both 2-D and 3-D SSB models. 
Mode skills are evaluated in terms of the explained variance 
We compared two 2-D SSB models (2D_US and 2D_MS) We 
fit calculte the explained variance withthe global data and 
resa indicate that: 2D_MS i better on a global seale. We then 
binned the crossover data into 10" atitude bands and computed 
the explained variance in each band. Results show that 2D. MS 
is better than 2D. US at most latitudes, We compared the 3-D 
‘SSB model (3D_USM) with the vo 2-D SSB models as well as 
the SSB correction in Jason-2 GDR product. Results indicate 
that the variance explained by 3D_USM increases by 1 32cm 
and 0.84 cn relative to that explained by 2D. US and 2D_MS 
respectively. Spatial evaluation of improvement shows that 


3D_USM is bener than both 2D US and 2D MS at al 
aides 

SLA is an imponan waditional parameter used im 
oceanographic investigations. We compared the Jaon-? derived 
SLA withthe tide gauge derived SLA at Corsica absolute 
‘calibration site- Results show thatthe variance of Jason-2 SLA 
bias from the tide gauges can be reduced by 0162-067 em 
When replacing the 2-D SSB estimates with the 3-D SSB 

he SSB correction is very important for the coastal 
altimetry. We tested the performance of the 3-D SSB estimates 
Jn the Coastal regions. The variance of the crossover SSH 
diterences as well as the variance of the along tack SSH 
‘difference was used to assess the 2-D and 3-D SSB estimates. 
Compared with the 2-D SSB estimates, the -D SSB estimates 
perform beter in the whole coastal regions. 

Unlike Travs 3-D SSB model, we used the mean wave 
period from ECMWF reanalysis project ERA-Imerim instead 
of that rom NOAA's WaveWatc model, We deseribed some 
agreements and disagreements between the two sources of the 
mean wave period. We nt compared the rwo 3-D SSB models 
using MWP from ERA-Imerim and Tm estimated from 
aveWatchà model respectively. Results show that the 3-D 
SSB model using MWP fom ERA-nirim has higher 
reduction in the total variance of the SSH differences. We then 
developed a -D SSB model using the peak wave period from 
WaveWatch3 venion 222 hindeast reanalysis product 
produced by NOAA/NCEPINWS and compared it with the 3-D 
SSB model using MWP fiom ERA-Interim. The 3-D SSB 
model using MWP frons ERA-Interim performs better than the 
3D SSB model using the peak wave period in the whole and at 
alt iud 

"The operational SSB models applied to the SSB correction 
{or Jason- 1 and Jason-2 are developed from SSH differences. 
Besides, the ECMWF products are widely used in altimeter 
level 2 processing. Both the dry oposphere path delay 
correction and the layered barometer correction use the daa 
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fiom ECMWF atmospheric pressure products. The ECMWF 
america weather prediction models may also provide a wet 
"ropospheri correction for Jason-1 and Jason-2 GDR products, 
Therefore. the method presented in this paper to estimate the 
[SSB is very practical for the SSB correction in radar altimeter 
measurements of sea level. The remis in this study are 
‘bluned with Jason-2 GDR daa The same method need to be 
applied to other altimeter mission data to insure thar the method 
presented in this study îs applicable to all data sets. 
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A. Pseudocode for LSMR 
LSMR is used io solve the systems of linear equations 
Ax =b Ire system is inconsistent, it solves the least-squares 
problem minjo- Ax, -To this paper, we use the LSMR to solve 
equation (39). B,, pand AY-B,g, are corresponding 10 A 
xand b respectively. We only give the basie steps for the 
LSMR algorithm. If the readers want to know more about the 
SM algorithm, they can refer to the Herature [6] or the web 
‘te hup/veb stanford. edu/groupSOLsoftwatfsni 
1 (italize) 

Bu, 
Beat asl d. 
2 For k=1.2.3..., repeat steps 3-6 
3. (Continue the bidiagonalization) 

B.u, SAV, =m, 
4. (Construct and apply rotation P, ) 


A'u, 


6. lUptaeh i, v) 
B an, p, Mpa 


v ott inp e 


ameta RR 
7. Compute s [ATE Js. estime n. cond) and 
ted if any o e sopping reri is sa 


B. The steps for deriving the SSB values with a FD 
nonparametric model 

1. Randomly extract N crossover points each time. In this paper, 
we choose N-80. 

2. Impose etx], With any reasonable value afg. we 
choose gp, =-12em in this paper. 

3. Compute gin (36) using the LSMR algorithm. 

4. Compute ex) over a regular grid of (U, SWH, MWP) by 


using 33). 
5. Repeat tps 1-4 for M mes In this paper, we choose MeS0, 


6. Compute x) the average of pix 
T. Subiact GU =0,SWH =O. MWF 
pre 


9s) from Bex) ata 


Remtaences 


LI) £ San and V. Zieni, Seu Sa a Rabar ety 
ta) ge any Z Zande, in be 
in ian sh Rac a WR 
(0) FP De Root: Nk Be oan Fes 
1) BS Vals AS Samm D An dE A an 
NS) AUN Seine Os te oa Din of Sree Eat and 
Sls in Soc Ra Soh spc o Rl 
iy ao Gye! Roc re VOTE 
I) Rage Y. 3. m and) N Marin "Te Ei of Se Wane 
Secs Ent at ped 
my Beha geile Vaca ud B Cha, Wal 
seni ryan econ maf Cope 
i| Tib DR Tumor B Capt ad D Vino 
“epee eret m ans rri Cd 
IT Rida D Tampa. B Chop and D. Vaden, 
eene ainiin Da any Eram of dpi 
ELLE pl kei d MR 
(10) W Aaina, K Wom, A" at 
re ice agi fun of Gye Rah De 
(un) PW Min K Wari, 1 Nae ad D- V. Arai, Pic 
‘tan eran Sh wee 
(us) Pea PT ie A Me Ca Sao rae 
iai ETE Nude DV. Hamat. D E ln nd UE Rey, 


201703.00307v1 


chinaXiv: 


IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING E 


ei 


m 


em 


ots ae ier ge 
M e ETT 
eae D erede peret 
FW Milla, D. V. Arii, P. Gaspar, KF Warnick; and 1 Soh 

aeri pone Juna of ipla Renard Dc, v. 
ec eer 
[eor it i 
Yann, N Tar BB Beckley, B. Se, FC 
ein bet ere 
Lee er TA ee 
3 Manga Lg Wine Means jir rbd 
S enia M Ya. Wn Theory and Cát Pr aper Being. 
Ifa and Gel oral oeil Mating an Ande 
Renner Dic ue AREE E 


(25) D.C L Fong ad M. Sanden, “Lane An bie Agi ir 
Bip ES M Uppal A 1. Semen, P ron, Pa, $ 
Toby, i aly "e ERA dim esse faim ad 
peace oh i nm yom” Quan rel f e 
A‘, Rule N ca and E Moss, ce Chl 
[D 

Cata cs c ois Matar Grey 3 P 


Mac ang sie he BS deg ia tie 


ETIN TOME 
pA 


Yanni à igi 
Technology. Nang Chin ad ie MS igi 
stl orp i an ae Ya 
sy un Dom Unely of Chin, Qin 

‘Buoy eiie m ctu indc. 


